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The nine Reactpr Stat ics Modules are designed to 
introduce students €b thq use of numerical methods and digital 
computers for calculation of neutron^flux distributions' in space and 
energy which are needed to calculate. criticality, power distribution, 
and fuel burn-up^ for both slow neutrtfn and fast neutron fission 
rsactors. The diffusion approximation is used for .the calculation of 
neutron transport. Collision probabilities are used td calculate the 
effect of heterogeneous lattices of ftiel and moderator , on/resonance 
absorption and slow neutron disadvantage faq&ors. In" this module, the 
three major problems of neutron diffusion in slow neutron Mission 
reactors are "reviewed, an energy-dependent diffusion equation, steady 
*state eoithermal multi-group equations, and a slow neutron diffusion 
equation' are introduced. The structure and organization of the 
remaining eight modules are outlined. These modules are intended to/ 
supplement textbooks and other lecture material generally available 
to students in their course work. It is assumed^ that students a^p4 " 
familiar with elementary nuclear structure, neutron-nuclei 
interactions, and introductory material on fission chain reactorsw 
(Author) 



* * * * * * * * * * * * * * * * * * * * in * * * * * * * «»t4K Jle^t 4t «|i 4|r4|r #t #r #r^t #r#i #r«|r #e 4r ilr^r^r^t «H 4> ^^t-^r «|ri|c % s|r ^tsli^ #r <K4i4ca|Ka 

* Reproductions supplied by EDBS are the best that can be made 

* from the oriainal document* 

** * * He * ** i|i * * * * * * * * ** * * * i|i i|t * * * * 4c4r4t4r4|e4B4e4^|i4r4i4t4i4K4K«|r4i«|c«|r4i4Ki|i «|r4e4i«|r 4e#c4>a|i 4b s|r 4t4fi|is|r 41 ^t^r^e^cs 

/ 



■ V 



INTRODUCTION TO REACTOR STATICS MODULES, ' RS-1 




Milton C. /Edlund 



4^ 



The University gratefully acknowledges the support of the. 
Division of Higher Education of the National Science Foundation 
for .support of this work performed under Grant GZ^2888 and the 
f support of Duke Power Company', North Carolina Power and Light Company, 

and Virginia Electric and Power Company. 



Project Director: Milton C. Edlund 



TABLE OF CONTENTS 
> - 

, Pa B e 

» 

1.1 The Three Major Problems of Neutron Diffusion In 

Slow Neutron Fission Reactors * 3 

1.2 Ar> Energy Dependent Diffusion Equation. • -\ * 5 

1.3 Steady State Epithermal)«Multi-Group Equations\ . .' • . 9 

1.4 Slow Neutron Diffusion Equation'. .< 16 

1.5 'Structure and Orgpnization'of the Reactor Statica Modules 17 



4*' 



< 



f 



INTRODUCTION TO REACTOR STATICS MODULES, RS-1 

. . . \ 

Preface 

The development of reactor physics shortly* after the discovery of nuclear 
fission in January 1939 is truly an elegant application of nuclear physics, and 
applied^ mathematics . The notion of a nuclear chain reaction wa$ intensively 



pursued 



after the disdovery of fission.^ Nuclear physicists at that^ time 
realize^ that neutrons would be emitted in the fission process.^. And since 
neutrons have no. electrical charge, they would easily penetrate the electric 
field surrounding a nucleus to interact with it and thereby serve. as the, 
carrier for a chain reaction. - 

During the early months of 1940 it; was s.hown that fission of uranium was 
primarily induced by slow.neutron absorption in the less abundant uranium -235 
isotope. With remarkable rapidity, the physical problems* of constructing a 
self-sustaining chain reaction in^rtatural uranium were identified and solved 
leading to the first nuclear fission chain reaction' at the University of 
Chicago on December 12, 1942. ^ 

2 ' 

Reactor physics developed during* the Manhattan Project as an inspired 

combination of some elements of mathematical physics and the definition and 

•^measurement <*f certain nuclear cross-sections and integral quantities related* 

to the mathematical models used in the design of the plutonium production 

reactors built at Hanford, Washington during World -War II. The models and 

analytical methods, which can be identified as* "Integral Reactor Theory," were" 

devised to minimize computations which were performed on desk calculators^ 

3 4 5 6 

Integral reactor theory is described in several textbooks, * ' ' which have 



5 



been used In nuclear engineering education during the last two decades. 

Two developments have led to a significant advance in the methods of 
reactor physics analysis. The fl^rst is the more accurate measurement of a 
large number of nuclear cross sections. The second is the development of high 
apeJd and large capacity digital computers. The latter permits the use of more 
sophisticated mathematical models for computation; the former provides the data 
required to make the more sophisticated models useful in the design of nuclear 
reactors. The use of digital computers also facilitates the exploration^of 
many core design variables and permits the solution of many problems which 
simply could not be- solved otherwise. 

> 

The Reactor ,Statics Modules are designed to introduce students to the use 
of .numerical methods and digital computers for calculation of neutron flux 
distributions in space and energy which are needed to calculate criticality, 
power distribution and fuel burnup for both slow neutron and fast neutron,, 
fission reactors. The diffusion approximation is used for the calculation of 
neutron transport. Collision probabilities are used to calculate the effect 
of heterogeneous lattices of fuel and moderator on resonance absorption and 
slow neutron disadvantage factors., V ^ 

These Modules are intended to supplement textbooks and other lecture 
material generally available to students in their course work. It is assumed, 
that students are familiar with elementary nuclear structure, neutron-nuclei 
interactions and introductory material on fission chain reactors at a level of 
Nuclear Reactor Engineering , S. Glasstone and A. Sesonske, Van Nostrand 
Reinhold Company, New York, N.Y. (1967). 
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i 1.1 The Three Major Problems of Neutron Diffusion in Slow ( . 
* Neutron Fission Reactors 

!• *>' - ' ' " 

There are three problems in diffusion theory, the solutions of whiph are 

s. ■ ■ • , ' ■ 

'necessary for the calculation pf the multiplication factor in slow. neutron 



reactors. The^ are. 



1) the slowing dowrf of fission heutrons to thermal energies' 

v 



without capture, 

2) thermalization "arJd diffusion of slow neutrons and 

3) the slowing down of fission neutrons taking into account 

^ t ^ e possibility of their capture in the resonance region. 

I ' ' ' . 

*v* fhe solution bf these problems ^provides the neutron distribution in energy 

s . and time throughout the reactor.* The neutron distribution and the macroscopic 

- p . r *■ / 

*£' cross sections give neutron-nuclei reaction rates, which are needed to calcu- 
late the multiplication factm: and power distributions. The neutron distribution 
is the primary dependent variable required fpr tl^e nuclear design of reactors. 

The slowing down of neutrons is due almost entirely to elastic scattering 
with moderator nuclei. In Tight nucMtei tfie first excited states lie about 1 Mev 
a^03(e /the ground state, and hence, .in£^a©1§Lc collisions are relatively unim- 

portant in the mod^ratfion process. -?articijlarly , since nuclei which are good 

\ ^ ft 

moderators by virtue of the elastic scawrering mechanism are of necessity light 

Elements* Inelastic scattering in the Heavy elements can reduce the energy of 

fast neuktons Gnly downj to about, Mev, which is the excitation energy of the 

first excited, state in heavy nuclei. This^process is important, however, in 

v 

determining the slowing down length in olcigue-packed uranium oxide fueled light 
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water reactors.. It is less important in gas -cooled graphite moderated reactors 

\ . ' ■ • ; . 

and is a primary process in fast breeder reactors. r t , 

Every elastic collisioh of a neutron with a ipoderator nucleus will reduce 
the N energy of. the neutrons until the kineticv^nergy of the ^eutron approaches 
an energy on the order of the thermal energy of the moderator atoms. At this 
Ctage, energy may be transferred from the moderator atoms to the neutron and 
the neutron energy will approach '"thermal energy 11 asymptotically. 

Above about 1 ev the moderator atoms can be considered to, be at rest 
compared to the motion of the neutrons and the calculation of the energy loss 
per scattering collision is an elementary problem in particle mechanics. 
Slowing' down above 1 ev is discussed, in Module RS-1. 

♦The last phase of the slowing down process, however, is quite difficult 
ttO analyze. Strictly speaking, these are no thermal neutrons in the sense of 
neZitrons being in thermal equilibrium with the moderator. Because' of the 1/v 
law of neutron absorption, the* absorption probability increases, at the end of 
the slowing down process 4nd neutrons are actually absorbed before attaining 
the energy distribution which corresponds to thermal equilibrium., The diffusion 
and energy distribution 'of slow neutrons (below -1 ev) is further complicated 
by the effect of chemical binding and, in, the case of crystal lattices such as n 
graphite, the effect of coherent scattering of low energy neutrons by many 
nuclei in the crystal. This problem is discussed in Module RS-5. 

The calculation of the capture of neutrons in the resonances of nuclei is 
complicated by the rapid variation of cross sections with energy. Nevertheless ,- 
the theory of this process is quite satisfactory and calculations of resonance 
capture dgn be made accurately and quickly with digital computers. This problem 

' j 

is discussed in Module RS-4. '3 



1«2 An Energy Dependent Diffusion Equation 

» 

Diffusion theory 18 widely used to calculate the neutron distributions 
needed for reactor design. The basic approximation of diffusion theory is 
Pick 1 8 law, which states that the neutron current per unit energy is proportional 
to the gradient of the neutron flux per unit energy; i.e., 

J(r, E, t) - - D(r, E) 7 *(r, E, t). (1.2.1) 

The diffusion coefficient is approximately given b^, 

D(r, E) = (1.2.2) 

Z (r, E) -) u (r, E) E (r , E) 
t ~~ L — i I — si ~~ 

■ i 

Z is the total macroscopic? cross section, £ . is the scattering cross .section 
t sx 

and is the average cosine of the scattering angle for the i th isotope. 
The summation is over all isotopes. 

Both the current and the flux are distribution functions in space, and 
energy, E. That is, the neutron current for neutrons in' the volume element d_r 
at r and having kinetic energies in the range E' to E + dE at time, t, is, 

r * . 

J(r, E, t)dr dE. 1 

Similarly, <f>(£> E, t) d£ dE gives the number of neutrons in djr at r_ with energies 
in the range E x to E + dE at time, t, multiplied by their speed, v = /2 E/m. 

Fick ! s law and the formula for the diffusion coefficient can be derived 

^ 7 

as a first order approximation from the exact neutron transport theory. The 
approximation is sufficiently accurate ffyr^ those reactor design problems in 
which the variation of the flux with r is small compared to a neutron mean free 
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path aiul the variation of flux with energy at each poultion la amall over the 
range of the average neutron energy change reuniting from a scattering collision. 

The diffuaion approximation can be rationalized by noting that the m imp lout 
relation between the vector field, J[, and the scalar field, which ia in- 
variant to rotation anc^ translation of the coordinate system, is precisely 
Equation (1.2.1). After all, the neutrons in a reactor don't "know" what 
coordinate system is being used by the analyst! 

Diffusion theory, however, is not sufficiently accurate for design re- 
quirements in the calculation of neutron distributions within a single lattice 
of fuel, clad and moderator. These distributions require more accurate neutron 
transport models. On the other hand, the migration of neutrons over a number 
of lattices can be treated by diffusion theory. The classical approach is to* 
divide the neutron distributions into "macroscopic 11 and "microscopic" components 
as illustrated in Figure 1.2.1. 

< 

Figure 1.2.1 



Fuel 



Microscopic Slovr Neutron Flux 
in Lattice Cell 




Macroscopic Slow Neutron Flux 
over Several Lattice Cells 
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Fuel Rods 



The microscopic dintrlbut Loiui , those within a fuel lattice coll, are computed 

using tho more accurate transport models to obtain coll average macroscopic 

cross sections. These crotia aoctlona are then uaed [\\ diffusion equations' to 

calculate the average neutron distributionra^over a numbcrt: of lattice^; I.e., 

the "macroscopic 11 distributions. An Introduction to the calculation of the 

"microscopic 11 distributions is presented in Modules RS-4 and 6. 

Returning to the "macroscopic" diffusion theory, the basic principle is 

v 

simply neutron conservation in an element djr d^at £ and E. To simplify the 
presentation, consider a single fissionable isotope such as uranium -235 mixed 
with a moderator consisting of a single isotope like graphite . Further, assume 
a dilute system in which the critical concentration o^ uranium -235 is so small 
that scattering by the uranium can be neglected. The most important neutron- 
nuclei interactions are fissions below about l'ev apd elastic scattering of 
neutrons by moderator nuclei. The time rate of change of the number of neutrons 
in dr dE is equal to the sum of x. 
1) the rate at which neutrons are scattered into dE by 
moderator nuclei 

and 



minus 



and 



2) the rate of production of fission neutrons 

3) the rate of leakage by diffusion 

4) the rate of removal by all interactions with the 
uranium and moderator. 



The number officii! I :«r lug noil Uilonu lu ilr illi 1 at lime, L , lu 

(y.i^)'* (iiK' .^) % * 

E fl (r jK 1 ) tfi the maorot&oplc ticutturlftH c^omii noctLon of moderator. Lot ( 

'a 

l»(IS f -►E)(lk 

be the probability that neutrons scattered with*' initial energy, K'.Yave 

' / 

final, energy in the range K to E + dE. The total number ot : neutronti scattered 

v \ 
into dE is obtained by integration with E' from 0 to 10 Mev and gives the first 

\ * . ; , 

term in the neutron balance, which is 

10 Mev 

dr dE j ^ 3 ( f I^ , ) * (r,E',t) P (E 1 + E)dE'. (1.2 J) 

- 0 

The production rate of fission neutrons is 

10 Mev ' 
. dr dE vX(E) j ^fCr^) <J> (r,E\t). (1.2.4) 
0 

X(E)dE is the fraction of fission neutrons emitted in the range f to E + dE, 
v is the number of fission neutrons produced per fission and E f (r,E') is the v - 
macroscopic fission cross section. 
The rate of leakage is 

- dr dE 7 • D(r,E) V<j> (r,E,t) (1-2.5) 
and the rate at which neutrons are removed by all interactions is , 

dr dE E t (r,E) <f> (r,E,t), (1.2.6) 
where E..(r,E) is yrfie total macroscopic cross section. 
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The time rate of change of the number of neutrons in dr. dE is ^ 



1 3*(r,E,t) 



- - ., dr - d % ? — — • ' (1 - 2 - 7) 

Neutron Conservation leajds to an energy dependent diffusion equation for 

the neutron flux per unit energy, which is 

; ,. •. * KMflev 

i a*(r, 



v 



f r,E,t) i C ' 

- V ■ D(r,E)7*(r,E l t)-E f (r,E)*(r,E,t)+ E ( r , E ' ) <f> (x, E \ t) P (E ' +E) dE 1 

dt L J S ' ' 



10 Mev 



+ vX(E) I I (r,E , )(}>(r,E , ,t)dE l . (1.2.7) 

J 0 



This is the basic form of the diffusion equation used in the development of 
the Reactor Statics and Reactor Dynamics Modules. Scattering by more than one 
isotope merely gives an additional term having the same form as the in-scattering 
term in Equation (1.2.7) for each isotope. Similarly, fission of more than one 
isotope yields additional fission neutron production terms. 

1.3 Steady State Epithermal Multi-Group Equations 

In the epithermal energy region, above ~1 ev, the energy exchange With 

moderator can be calculated as though the moderator nuclei are at rest relative 

Ato neutrons. Also, the elastic scattering with low mass moderators is isotropic 
fo*l energies up to the Mev region. The slight anistropy at higher energies will 

be neglected. From elementary mechanics (cf., Glasstone and Sesonske, pp. 127- 

132), the probability that neutrons scattered with initial energy, E f , have 
final energy in the range E to E + dE is 

^ P(EU- E) dE - dE , gt for E < E' < - (1.3.1) 

* (l-a)E - a 

= 0 otherwise j?'?*^, 



1 O 

X KJ 
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jwtiere 



A - 1 




and „ / 

Mass of scattering nucleus? 

^ C3 ■ ■ M - ■ t ■ . # 

Mass of neutron 

For example, if the moderator is carbon, a - 0.716 and neutrons can be scattered 
into dE a«t E only if their initial energies are in the range of E to 1.40 E. 
Note, also, that Equation (1.3.1) states that neutrons are uniformly scattered 
over the range E' to aE'. [ 

The epithermal energy region can be divided into a number of contiguous 
energy groups. The g th group has the energy range, E - E let g - 1 denote 
highest energy group, g = 2 the next highest etc. Generally, a reactor core and 
reflector or blanket can be divided into a relatively small number of regions, 
within each of which, the macroscopic cross^ sections are independent of jr. Thus, 
setting 3<t>/Bt - 0, using P(E f -+ E) given by Equation (1.3.1) and integrating 
Equation (1.2.7) with respect to E over the range, E - E , gives for each such 

reactor region, 
E E E^ E/a 

Jg? r R r g r dE* 

D(E) V% (r,E)dE - I E fc (E) * (r,E)dE + / dE I L g (E 1 ) * <r,E f ) (1 _ a) ^ 

E s-i . 

E 10 Mev 

g 



X(E)dE J 



+ v | X(E)dE J E f (E') $ (jr.E^dE' = 0. (1.3.2) 



E , ~ 0 

Next, assume that $(jr t E) is separable over the range E g - E g _ x throughout the 
region; i.e., let 

1 * 

<Kjr,E) - G(R) F(E) . 



Then 



t if 
E « r > E 

g • ' ■ g 



f D(E)V 2 <(>^T,E)dE - V 2 q(r) / J D(E)F(E)dE. (1.3.3) 



E , / - E 

g-1 / g-1 

(g) 

Defining D as the spectrum average value of D(E).; I.e., 



E 

g 'V. 



f D(E)F(E)dE 

J 1 



D (6) = , ' (1.3.4) 

" g 

F(E>dE 



/ 



g-1 



and noting that the group flux is 
v 



E g E g 

* g (r) = J <t>(r,E>dE - G(r) J F(E)dE, (1.3.5) 

E E 
g-1- g-1 



the leakage term of the group equation becomes 

E 



Similarly, 



^ g V D(E)V 2 <t,(r,E)dE - D (g) V 2 4> (r) . (1.3.6) 



E 

g 



/ E t (EH(r, 



E)dE = z/ 8) * (r) , (1.3.7) 



E C * - 

g-1 
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with 



■r. 



(E)F(E)dE 

E ( 8 } - ^ . , - (1.3.8) 



Also, 




E E/a g . 



,f dE f E(E') » (r.E') 3 d *' t ■ IlM^^d). (1.3.9) 

E r (j * 8) is the transfer cross section fromgroup j to group g f which when 
multiplied by * (r) gives the rate at whi^^putrons are -scattered from group 

j to group g per unit volume at £. The tranafer cross section is , 

<* 

M dEp E 8 <E')F(E') TI ^ ¥7 
J E . J E . 

I 

/ 



Er (j . g) . Jj^L ^fcl (1.3.10) 

g 

F(E)dE 



Vi 

In the transfer cross section, the upper limit of the integration over E 1 is 

E /a, since P(E f E) vanishes for higher energies, 
g 

If inelastic scattering is important, Equation (1.2.7) will contain a term 
of the form, 



10 Mev 



/ 



Z, (E f )F(E')K(E' + E)dE', 
in 1 p 

E ^ ^ 



- 13 - V * : ^ 
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where K(E f + E) is the transfer probability for, inelastic scattering and - 

L. (E f ) is the totkl inelastic cross section f6r each inelastic scattferer. The 
in 

inelastic transfer cross section is calculated in the same way« as the elastic 
transfer cross section. 



The fission neutron production term in Equation (L. 3.2) becomes 

N \ 



vX (g) £ £ (g) # (r)< 



The summation is over all groups including a thermal or slow neutron gr^oup, 

The £^ g ^ is the' average value of over the spectrum F(E) and ( , ' 

E „ ' 

X^^o I X(E)dE. 



E . 

1 g-1 
The epithermal multi-group equations, thus have the form, » 

8 N " 

D (g) V* g (r)-Z t (g) $(r) + X! V^ 8) *j (l)^ 8 ^ z f (8> * g <£> =0 - (1-5.12) 

J-l g»l " . 

• g = 1, 2 . . . N - 1. 

N is the thermal or slow neutron group which is discussed in section 1.4. 

Turning, now, to the calculation of the epithermal group constants, the 

i 

simplest approximation is to use an infinite medium spectrum. Equation (1.2.7) 

? 

reduces to 

E/a 

- Z t (E)F(E) + j £ g (E , )F(E , ) (1 ? E ^i + X(E) - 0. (1.3.13) 
E 

F(E) is the neutron spectrum normalized to one fission neutron produced per 
unit volume and unit time; i.e., 



1 ^ 



10 Mev ' . r \ 

vj E £ (E , )#(E , )«' - 1. ~ cAl4) 

° y 
The accuracy of this approximation for the flux spectrum within each group m 

increases with decreasing energ^width of the group. Thus, for a" fairly large 

number of groups, say 20 to 30, F(E) can be obtained by numerical solution of 

Equation (1.3.13) . • • £ 

Another approximation which is widely used is based on asymptotic reactor 

theory. If the cross sections are independent of r over a region which is 

large comparejl^tfo^a group dif fusion '"length, defined by, ^ 




L (1 - 3 ' 15) 
g V v (6) 



the neutron flux energy distribution wjill assume an asymptotic distribution' 
which is Separable in spade and energy over a region which- is about L away 

s 

from ii$frer£a , c s es with other regions. , - ^ ' 7 



As before', <J>(r,E) - G(_r) F(E). Furthermore, G(R) will satisfy the Helmholtz 



equation, 



V 2 G(r) + B 2 G(r) « 0. (1.3.16) 



The B 2 is the material buckling (cf., Glasstone ^n^ Sesonske, pp. 157-162). 

It is a function : $f the macroscopic cross sections and multiplying and scattering 

properties of the medium. If the infinite medium multiplication facto^^^,^ is 

2 2 
greater than zero, B will be positive. Conversely, if k^ < 0, B will be 

\ 

negative. ? 



Asymptotic reactor theory can be used for the calculation o'f th$ critical 
sizfe of a tfare gebme^^cslly 'simple homogeneous reactor. The spatial distri- 
bution of the flux per unit energy for all neutron energies / satisfies the - 
Helmholtz equation, 



V 2 G(r) + B 2 G(r) * 0, 



(1.3.17) 



where B is the geometric buckling which is determined by setting G(_r) = 0 pnf 
the extrapolated surface of the reactor. For example, a spherical reactor o& 
extrapolated radius, R p would have the spatial distribution, 



G(r) 



sin r 



and 



Vi 



/ 



.* v 



If the material buckling is equal to the geometric buckling, the re^tjor ^is 
critical, . « 

Xhe asymptotic reactor theory approximation for computing the neutron;, 
spectrum yields ttfe following equation for F(E) , 

E/a 



-[D(E)B 2 +I t (E)] F(E) 



V E? > F < E '>(1^ 



T + X(E) = 0. (1.3.18) 



The group constants are computed using this solution for F(E). 

The choice of the buckling to use in Equation (1.3.18) is made by trial 
and error. A guess is made, the spectrum and multi-group constants calculated 
whi£h are then used to calculate the spatial distribution of the group fluxes. 



1 n 

jl yJ 



Inspection of the group flaxes will indicate whether or not the average leakage, 
-D^B 2 <t> , over the region is consistent with the value of B^. . If not, a • 
different value of if can be selected and the process repeated. 

This approximation is the diffusion theory analog of the consistent P-l 

*~ ' ^ 
8 

and B-l approximations used in reactor design codes . 

In the neutron slowing dpwn module, RS-3, resonance absorption Is treated 

-'-'( - ' K / 

separately as difecussed Tn Module RS-4. Neutron absorption is divi<l£d into a 

\ * 'M ;<r 

smooth component tyhich is averaged over the calculated A&tra. and a resonance 

component. Resonance escape probabilities are comp^i^06x . e^& group. < 



■'l.4 Slow Neutron Diffusion Equation 



* In the thermalization energy region beloV a^iit f ev energy exchanges 
between neutrons and moderator are quite complicated as discussed in section 

1.1. However, because of the strong enetgy coupling between the moderator and 

'la 

neutrons, the medium imposes a definite energy distribution on the neutrons 

„ h 9 

which can be calculated by several more or less sophisticated models . This 

• * * 

distribution is used to calculate the slow or "thermal" group constants for a 
ohe group diffusion equation. * ^ 

I The steady state conservation equation is ^ ^ 

• ) 

D V 2 * (r,t) - I * (£,t) + q(r,t) ■ 0. (1.4.1) 



mi 



The slow group flux is the integral of the distribution from zero the maximum 



energy of the thermalization region, E c ; i.e., 

E 



* (r,t) = I <t> (r,E,t)dE (1.4.2) 
s is 



o 



2° 



• i At* ' ■ ts 



^ q(r,t) is \*ie* source of nex^rbns^ slowing down beKw E c and is given by 



4 E /a 



E 



E - aE 1 

c 



(1 - a)E' 



dE f . (1.4.3) 



Thi's equation for q can Jje derived simply by inspection. Neutrons at E' are 
scattered uniformly over the energy interval, (1 - a)E' . The fraction that are 
scattered below E ' is (E c - ctE 1 )/(l - aE f ) . ^ 

The slow neutron diffusion coefficient and absorption cross section are 
average values weigfiteji by the slow neutron flux energy distribution and <D g 

: ■ 

i^^the total slow neutron flux, *"^> , * 

The calculation- of sUw neutron /spectra using the hydrogen gas model is 

v 4 

discussed iK Modules RS-5 and 6. 



1.5 Structure and Organization of the Reactor Statics Modules * 
The following Reactor Statics Modules are: 



RS Number Title 

2 One Group Neutron Diffusion Theory \ 

3 Neutron Slowing Down and Epithermal Group Constants 

4 Resonance Absorption 

5 Slow Neutron Disadvantage Factors 

6 Slow Neutron Spectra and Group Constants * 

7 Three Group Criticality Program 

8 Multi-Group Constants for Fast Reactors 
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RS-2 introduces the student to the concept of the critical eigenvalue, k, 
the discretization of a one group diffusion equation and the numerical solution 
of the one group reactor mbdel for k and the. spatial neutron flux distribution. 

For slow neutron reactor calculations, the neutron energy range is divided 

into three groups, fast, resonance and slow. The division follows naturally 

from the different neutron-nuclfei interactions in each energy group. Thus, all 

fission neutrons are produced with energies in the fast group. Also, neutrons 

cannot be removed from the fast group by inelastic scattering. The resonance 

group contains all resolved resonances and is sufficiently wide so that neutron 

transfer by elastic scattei4ng from the fast group to the slow group can be 

neglected. Finally, the entire thermalization process is confined to the slow 

group. Their energy "rangd and lethargy widths are given in Table 1.5.1. They 

are called macro-groups to distinguish them from the micro-groups used in 

f ^ 
neutron spectrum calculations. 

Table 1.5.1 
Slow Neutron Reactor Macro-Groups 
Macro-Group Energy Range (ev) Lethargy Width 

1 7 

Fast 5.531 x 10 J - 1.000 x 10 7.50 



Resonance 



0.640 - 5.531 x 10 3 9.06 



Slow 0.000 - 0.64 

The group constants can be calculated for each reactor region using RS-3 
and 4 for the fast and resonance group and RS-5 and 6 for the slow group. One 
dimensional criticality and neutron distribution problems can be solved using 



- 19 - 



these group constants as input into the three -group criticality program 
presented in RS-8. 

The last module includes a separate .program for the calculations of 
neutron spectra and group constants for use in fast breeder reactor calculations. 
It can be used with the criticality programs in RS-2 and RS-8, 

A typical ^criticality and power distribution calculation for a low enriched 
uranium oxide slow neutron reactor using the three group model can be made using 
the following procedure, Fi^rst, the reactor is divided into regions, within 

i i 

which material properties are constant, for each region, the isotopic number 
densities in the fuel, cJTad ancf moderator are calculated. The fuel lattice cell 
geometry and the average temperature of* the fuel is calculated' using the 

■ - - 7 ' 

Thermal-Hydraulics Module, TH-1. This provides input data for each reactor 



region for the following computer programs; 

\ * 

Module Program Name Output 

RS-3 FARCON Fast and Resonance Group Constants 

RS-6 SLOWCON Slow Group Constants 

This output information is used as input data to the three group-one 
dimensional code, ODMUG, presented in RS-7. ODMUG can be used to compute the 
^critical eigenvalue, k, and the uniform slow neutron poison cross section 
required to make the reactor critical. The output of ODMUG includes; 

1) Critical eigenvalue of "unpoisoned" reactor, 

2) Poison cross section needed for criticality, 

3) Spatial flux distribution, 



n *7 
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4) Spatial power distribution and 

5) Total power generated and power density in each region. 
The sequence of calculations is shown in the following diagram. 



Physical Input 

U0^ fuel rod radius 

Fuel clad thickness 

Fuel lattice cell dimension 

Average fuel temperature 

Isotopic number densities 

Region size and geometry 



FARCON 



SL0WC0N 



ODMUG 



THREE-GROUP CONSTANTS 

— t 



V 



V 



Output 



Fast reactor calculations can be made by using FASTCON to generate group 
constants for either ODOG or ODMUG. 
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